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ABSTRACT 

We present here new results on the space density of rich, optically-selected, clusters of galaxies at 
low redshift [z < 0.15). These results are based on the application of the matched filter cluster- 
finding algorithm (as outlined by Postman et al. 1996 and Kawasaki et al. 1998) to 1067deg^ of 
the Edinburgh/Durham Southern Galaxy Catalogue (EDSGC). This is the first major application of 
this methodology at low redshift and in total, we have detected 2109 clusters above a richness cut- 
off of Rjn > 50 (or Ac; > 10; Postman et al. 1996). This new catalogue of clusters is known as the 
Edinburgh/Durham Cluster Catalogue II (or EDCCII). We have used extensive Monte Carlo simulations 
to define the detection thresholds for our algorithm, to measure the effective area of the EDCCII and 
to determine our spurious detection rate. These simulations have shown that our detection efficiency 
is strongly correlated with the presence of large-scale structure in the EDSGC data. We believe this 
is due to the assumption of a flat, uniform background in the matched filter algorithm. Using these 
simulations, we are able to compute the space density of clusters in this new survey. We find 83.5t.lg%^ x 
lO^^h-^Mpc-^ for 100 < Rrn < 200 (A^z ~ 20) systems, lO.l^^g^ x 10~^h~^Mpc~^ for 200 < i?„, < 

400 (Ac; ~ 40) systems and 2.3t2.3 x 10"^ ^"^ Mpc'^ for i?„ > 400 (A^; > 80) systems. These three 
richness bands roughly correspond to Abell Richness Classes 0, 1 and > 2 respectfully. These new 
measurements of the local space density of clusters are in agreement with those found at higher redshift 
(0.2 < Zest < 0.6) in the Palomar Distant Cluster Survey (PDCS; Postman et al. 1996 & Holden et 
al. 1999) and therefore, removes one of the major uncertainties associated with the PDCS as it had 
previously detected a factor of 5 ± 2 more clusters at high redshift than expected compared to the space 
density of low redshift Abell clusters. This discrepancy is now lessened and, at worst, is only a factor 
of 4^4^. This result illustrates the need to use the same cluster-finding algorithm at both high and 
low redshift to avoid such apparent discrepancies. We also confirm that the space density of clusters 
remains nearly constant out to z ~ 0.6 in agreement with previous optical and X-ray measurements 
of the space density of clusters (Couch et al. 1991; Postman et al. 1996; Ebeling et al. 1997; Nichol 
et al. 1999). Finally, we have compared the EDCCII with the Abell catalogue. We detect nearly 60% 
of all Abell clusters in the EDCCII area regardless of their Abell Richness and Distance Classes. For 
clusters in common between the two surveys, we find no strong correlation between the two richness 
estimates in agreement with the work of Lumsden et al. (1992). In comparison, ~ 90% of the EDCCII 
systems are new, although a majority of them have a richness lower than an Abell Richness Class of 
and therefore, would be below Abell's original selection criteria. However, we do detect 143 new clusters 
with Rm > 100 (which corresponds to a Richness Class of greater than, or equal to, 0) that are not in 
the Abell catalogue i.e. 63% of the rich EDCCII systems. These numbers lend credence to the idea that 
the Abell catalogue may be incomplete, especially at lower richnesses. 

Subject headings: cosmology: observations — galaxies: clusters: general — galaxies: evolution - 
surveys 



1. INTRODUCTION 

Clusters of galaxies play a key role in tracing the dis- 
tribution and evolution of mass in the universe (see, for 
example, Guzzo et al. 1992; Postman et al. 1992; Nichol 
et al. 1992, Dalton et al. 1992; Bahcall & Soneria 1983; 
Reichart et al. 1999). Until recently, such studies have 
been based on catalogues of clusters constructed from vi- 
sual scans of photographic plates e.g. the Abell catalogue 
(Abell 1958; Gunn, Hoessel & Oke 1986; AbeU et al. 1989; 
Couch et al. 1991). However, during the past decade, 
there has been considerable progress in the construction 
of automated catalogues of clusters and groups that pos- 
sess objective selection criteria. Such work includes clus- 
ter catalogues selected from digitized photographic mate- 



rial (Dodd & MacGillivray 1986; Lumsden et al. 1992; 
Dalton et al. 1994), from X-ray surveys (Kowalski et al. 
1984; Ebeling et al. 1997; de Grandi et al. 1999; Gioia et 
al. 1990; Nichol et al. 1997 & 1999; Rosati et al. 1998; 
Burke et al. 1997; Jones et al. 1998; Romer et al. 1999) 
and large-area optical CCD surveys (Postman et al. 1996; 
Lidman & Peterson 1996; Zaritsky et al. 1997; Olsen et 
al. 1999). 

There has also been great progress in the development 
of new cluster-finding algorithms. The first automated 
cluster catalogues used simple variants on the "peak- 
finding" algorithm (Lumsden et al. 1992) or the percola- 
tion method (Dalton et al. 1994). In recent years, several 
new, more sophisticated, algorithms have become avail- 
able including the matched filter algorithm - in several 
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Fig. 1. — The 10° x 10° test region outlined in Section ^ and used for our simulations. Cluster detections within 5 X 8c (at z = 0.05 this is 
0.35°) of the edges of this data were discarded. 



different flavors (Postman et al. f996; Kawasaki at al. 
1998; Kepner et al. 1999; Schuecker & Bohringer 1998) -, 
the wavelet-filter (Slezak et al. 1990), the "photometric 
redshift" method (Kodama et al. 1999), voronoi tessella- 
tions (Ramella et al. 1998) and the "density-morphology" 
relationship (Ostrander et al. 1998). The level of sophisti- 
cation of these algorithms has increased in anticipation of 
high quality CCD survey data e.g. the Sloan Digital Sky 
Survey (SDSS; Gunn et al. 1998), 

The early automated catalogues of optically-selected 
clusters have produced two important results. First, Post- 
man et al. (1996) and Lumsden et al. (1992) both find 
evidence for a higher space density of clusters than that 
seen in the Abell Catalogue. For example. Postman et al. 
(1996) finds that the measured space density of clusters in 
the Palomar Distant Cluster Survey (PDCS) is a factor of 

5 ± 2 greater than that implied from the Abell catalogue. 
Second, the space density of PDCS clusters remains con- 
stant between z = 0.2 and z — 0.6, in agreement with the 
earlier work of Couch et al. (1991) and has been confirmed 
recently by Holden et al. (1999). If true, these results 
can be used to place strong constraints on the underlying 
galaxy evolution model {e.g. CDM) and measurements of 
the cosmological parameters as and ilo (see Bahcall, Fan 

6 Cen 1997; Reichart et al. 1999; Holden et al. 1999). 
To solidify these initial results, larger catalogues of clus- 
ters are required. Moreover, it is becoming increasingly 
clear that we need to compare these different cluster cata- 
logues to help verify results and expand the redshift range 
over which we can study the cluster distribution. To date 
however, there has been little cross-comparison between 
these different cluster catalogues. Foremost, the relation- 
ship between the X-ray and optical catalogues of clusters 
remains unclear (see Holden et al. 1997; Briel & Henry 
1993; Bower et al. 1997). In the optical domain, different 
catalogues have used different cluster finding algorithms 



thus making it very difficult to cross-calibrate catalogues 
and methods and thus verify results. This is illustrated 
by the fact that although both Lumsden et al. (1992) and 
Postman et al. (1996) find a higher space density than the 
Abell catalogue, the PDCS finds 5 times as many clus- 
ters per unit volume, while the Edinburgh-Durham Clus- 
ter Catalogue (EDCC; Lumdsen et al. 1992) only finds 
twice as many clusters per unit volume as Abell. There- 
fore, it is impossible to fairly compare the EDCC and the 
PDCS even though they are both objective, automated 
catalogues of clusters. 

In this paper, we set out to rectify this problem by run- 
ning a variant of the PDCS cluster-finding algorithm on 
the same galaxy data as used by Lumsden et al. (1992) 
in the construction of the EDCC. The main aim of this 
project is to provide a coherent set of cluster data that 
spans from z ~ 0.05 - the lower redshift limit of the 
EDSGC - to z ~ 0.6 - the upper completeness limit of the 
PDCS. In addition to using a similar algorithm as Postman 
et al. (1996), we have performed a large number of Monte 
Carlo simulations to assess the completeness limit, and 
contamination rate, of this new EDCC cluster catalogue. 
This is the first major application of the matched filter 
algorithm to low redshift galaxy data, however, it is only 
the first of many such surveys presently underway e.g. the 
SDSS, DeepRange (Postman et al. 1998), DPOSS (Gal et 
al. 1999), COSMOS (Schuecker & Bochringer 1998) and 
the CCD survey of Zaritsky et al. (1997). 

In Section ^, we discuss the EDSGC catalogue and the 
matched filter detection algorithm. In Section ^, we out- 
line the methodology used to detect our cluster candidates 
and discuss in detail the Monte Carlo simulations we per- 
formed to determine our detection thresholds, the effective 
area of our new cluster survey and our spurious detection 
rate. Readers interested in just the results of this sur- 
vey may wish to concentrate on Section of this paper 
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which presents our space density results. In Section 5, we 
discuss these results in comparison with the PDCS and 
Abell cluster catalogues. Throughout this letter, we use 
Ho = 100 /ikms^^Mpc^^ and qo = 0.5 unless otherwise 
stated. 

2. THE EDINBURGH-DURHAM SOUTHERN GALAXY 
CATALOGUE 

The Edinburgh-Durham Southern Galaxy Catalogue 
(EDSGC) has been discussed in detail in Heydon- 
Dumbleton et al. (1989), Lumsden et al. (1992), Colhns 
et al. (1992), and Colhns, Nichol, and Lumsden (2000). 
However, for consistency, we include here a brief discuss 
of the salient points of this catalogue. 

The whole EDSGC comprises of 1.5 million galaxies 
brighter than bj — 21.5 covering an area of ~ 1100 squared 
degrees centered on the South Galactic Pole, spanning 
90 degrees in Right Ascension and 20 degrees in decli- 
nation. The catalogue is 95% complete to bj = 20 and 
has < 10% stellar contamination. The catalogue was con- 
structed from COSMOS (a microdensitometer) scans of 
UK Schimdt Illa-J photographic survey plates and as pho- 
tometrically calibrated using 30 CCD sequences positioned 
in a "checkerboard fashion". From the EDSGC, Lums- 
den et al. (1992) detected 733 galaxy overdensities using 
a simple "peak-finding" algorithm. Collins et al. (1995) 
presents 777 redshifts measurements within 94 clusters and 
this redshift sample has been used to study the large-scale 
distribution of nearby clusters (see Nichol et al. 1992, 
Guzzo et al. 1992 & Martin et al. 1995) as well as the 
cluster luminosity function (Lumsden et al. 1997). 

For the analysis discussion in Section |^, we restrict our- 
selves to a small 10° x 10° random subregion of the EDSGC 
centered at OO'"^'^ 30™™'' in Right Ascension and -33° in 
Declination. We also restricted the magnitude range to 
15 < bj < 20.5 to remain as complete as possible (see 
Collins, Nichol & Lumsden 2000). These cuts resulted in 
a total of 41171 galaxies which is significantly smaller than 
the whole EDSGC. This test data is shown in Figure |l|. For 
the space density results presented in Section ^, we used 
to all galaxies in the magnitude range 15 < bj < 20.5 
and in a coordinate range of 22^''*' < a < 3.3^''*' and 
-42° <S < -23° (1067deg2) which gave us 627260 galax- 
ies in total. 

We note here that the EDSGC has previously been used 
to construct an objective catalogue of clusters of galaxies 
(see Lumsden et al. 1992). However, in this prior analysis, 
only a simplistic "peak-finding" algorithm was used to find 
candidate systems for redshift follow-up. Given the recent 
advances in cluster-finding algorithms, we decided it was 
prudent to repeat the analysis which we discuss herein. 
We stress however that this does not undermine the scien- 
tific integrity of the original EDCC catalogue (Lumsden et 
al. 1992) and results derived from it (Nichol et al. 1992; 
Collins et al. 1995; Martin et al. 1995). In this present 
work, we simply wish to analyze the low redshift cluster 
population using the same techniques as presently used 
at high redshift {i.e. PDCS). For the sake of consistency, 
we call this new cluster catalogue the Edinburgh/Durham 
Cluster Catalogue II (EDCCII). 



3.1. The Matched Filter 

In recent years, there has been considerable progress in 
the development of new, automated cluster-finding algo- 
rithms (see Lumsden et al. 1992; Dalton et al. 1994; Post- 
man et al. 1996; Kawasaki et al. 1998; Kepner et al. 1999; 
Schuecker & Bohringer 1998; Slezak et al. 1990; Kodama 
et al. 1999; Ramella et al. 1998; Ostrander et al. 1998). 
In this paper, we focus our attention on the matched fil- 
ter algorithm since we are interested in directly comparing 
our results to those of Postman et al. (1996). 

We have based our matched filter algorithm on the pro- 
cedure outlined by Kawasaki et al. (1998) which compares 
the galaxy distribution around any point on the sky to a 
cluster model plus a background (sec Eqn. 1 in Kawasaki 
et al. 1998). The parameters of this cluster model are 
given in Eqns. 2, 3 & 4 of Kawasaki et al. (1998) and are 
the cluster surface density profile, the intrinsic richness of 
the cluster {N in Eqn 1. of Kawasaki et al. 1998), the clus- 
ter and field luminosity functions and the surface density 
of background galaxies (cr/). For the analysis discussed 
herein, we have modeled our clusters as a spherically- 
symmetrical, isothermal surface density profile {i.e. a King 
profile with rcore = 170 kpc and /? = |) combined with 
a Schechter luminosity function (with M^, — —20.12 and 

a = —1.25; see Lumsden et al. 1997). The background 
galaxy distribution is modeled as a flat surface density of 
galaxies of cr/ = 583775 galaxies per steradian which is 
the measured average surface density of galaxies in the 
EDSGC in the magnitude range 15 < bj < 20.5. For the 
field luminosity function {9 / (m) in Eqn. 3 of Kawasaki et 
al. 1998), we used a Schtecter function with Af^* = —19.5 
and an a — —1.1 (see Loveday et al. 1992). 

The physical model for the cluster plus background 
is converted to an observational model using the stan- 
dard cosmological redshift -distance relationships (and k- 
corrections) and therefore, the model is only a function 
of the intrinsic richness of the cluster and its redshift. 
This observed model is convolved with the data and a 
likelihood assigned for each point on the sky which is pro- 
portional to the quality of fit of the cluster model to the 
observed galaxy distribution given Possion statistics (see 
Eqns 6 & 7 of Kawasaki et al. 1998). One can then maxi- 
mized the likelihood by varying the two parameters of the 
model i.e. richness and redshift. Computationally, this 
was achieved by overlaying the EDSGC with a grid (pixel 
scale of Oc/3), and comparing each pixel, in this grid, with 
the matched filter model as a function of cluster redshift 
and richness {Rm)- For the results presented herein, we 
varied the matched filter redshift (zest) from 0.05 to 0.15 
{Szest = 0.0025) and used richness estimates of Rm = 50, 
100, 200 and 400. This resulted in an array of likelihood 
and richness maps (see Kawasaki et al. 1998) from which 
we must select our cluster candidates. 

3.2. Monte Carlo Simulations 

In this section, we outline our Monte Carlo simulations 
which were used to determine the cluster detection thresh- 
olds as well as to estimate our spurious cluster detection 
rate and overall cluster detection efficiency. 



3. METHODOLOGY 



3.2.1. Model of a Cluster 
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Table 1 

Richness threshold values (RT) determined for various redshifts and cluster richnesses 
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Table 2 

Log likelihood threshold values (LT) determined for various redshifts and cluster richnesses 



For our Monte Carlo simulations, we must first create 
an artificial cluster. We used a spherically-symmetrical, 
isothermal surface density profile with a cluster core ra- 
dius {tc) of 170 kpc, a cutoff radius of 5rc and a Schtecter 
luminosity function with a = —1.25 and Mg = —20.12 
(we restricted ourselves to an absolute magnitude range of 
M^^ ±5). This artificial cluster was then redshifted appro- 
priately using the standard cosmological relations and 

m = M + (42.384 + 5 log z) + (4.14z - 0.44z^), (1) 

to convert absolute magnitude to apparent bj magni- 
tude (see Lumsden et al. 1997). To match the magni- 
tude range covered by the EDSGC survey, we then re- 
moved all galaxies outside the magnitude range of 15 < 
bj < 20.5. Four different intrinsic richnesses of artifi- 
cial cluster were used in our Monte Carlo simulations; 
i?„ = 50, 100, 200, &400. Unless otherwise stated, each 
artificial cluster was unique, with the galaxies distributed 
at random according to the angular and luminosity distri- 
butions given above. 

We note here that no attempt was made to simulate the 
density-morphology relationship or to allocate a cD-type 
galaxy at the cluster core. Moreover, we did not change 
the shape or parameters of our artificial clusters during 
our simulations. Such simulations would have been com- 
putational intensive but are clearly needed in the future. 

more would have resulted in merging of the artificial 
clusters, while adding fewer clusters would have made the 
simulations laborious. We then varied the thresholds, both 
in richness (RT) and likelihood (LT), and computed for 
each combination the number of artificial clusters detected 
as well as the total number of clusters detected above these 
thresholds. 

It was then necessary to weight these detections to de- 
termine the optimal thresholds. Clearly, we wish to rule 
out obvious cases i.e. detecting 1 artificial cluster while 



detecting hundreds of other systems within the EDSGC 
(as most of these detections will be either lower richness 
clusters or spurious detections). Unfortunately, most cases 
we encountered were more subtle than this i.e. detecting 
10 artificial clusters within a total of 25 detections, or de- 
tecting 15 artificial clusters within a total 40 detections. 
To help differentiate between these intermediate cases, we 
employed an analytical method which we outline below. 

We defined the number of detections of our artificial 
clusters to be x and the number of total detections to be 
y. Both X and y are functions of RT and LT and increase 
with decreasing thresholds. We therefore defined a func- 
tion {T{x,y)) which satisfied the following boundary con- 
ditions; T(0, y)= 0, limy^oo T(x, y) = and the maximum 
of T{x, y) is at T(20, 20). The functional form of T(a;, y) is 
therefore important since it is now our weighting scheme. 
We estimated T(a;, y) by running our cluster-finding algo- 
rithm over many thousands of realizations of our simula- 
tions and varying RT and LT to create a Monte Carlo es- 
timate for T(x,y). From these data, we empirically found 

that T(a;, y)= was the best functional form for this 
data and adopted it as our weighting scheme. Using this 
relationship, we then determined the optimal RT and LT 
thresholds as a function of redshift and intrinsic richness. 
Our thresholds are given in Tables |l] and || (we used linear 
interpolation between these values if necessary). 

3.2.2. Locating Cluster Candidates 

In this subsection, we review our procedure for identi- 
fying a unique cluster candidate above the thresholds out- 
lined above. We first designate all pixels in our likelihood 
and richness maps that satisfy our thresholds, >RT and 
>LT, as "active" pixels. Obviously, a single cluster can- 
didate will create multiple active pixels so we must group 
these pixels together into single detection. This is achieved 
by searching for peaks in the distribution of active pixels 




Fig. 2. — The averaged measured detection efficiency for clusters of various richness as a function of redshift. Clockwise from the upper 
left-hand panel, the graphs represent clusters of richnesses Rm, = 400, 200, 100, and 50. 



i. e. we look for a active pixel whose height is greater than 
any other active pixel within a radius of 29c of the peak. 
We then group together all active points in this area to 
create a unique cluster candidate. If a single active pixel 
is a peak by default - i.e. there are no other active pixels 
within 20c of it - we disregard this peak and do not include 
it in our final analysis. Such isolated peaks are unlikely 
to be caused by real clusters. Finally, we compute the 
weighted mean of all active points grouped together as a 
single cluster to determine the most likely cluster candi- 
date centroid. 

3.2.3. Determining Detection Efficiencies 

In addition to determining our detection thresholds, our 
Monte Carlo simulations were used to estimate our detec- 
tion efficiency and the effective area of our cluster search. 
This is an important aspect of any cosmological survey as 
it allows us to determine the volume sampled by the ED- 
CCII and thus measure the space density of clusters from 
the catalogue. Previous applications of automated optical 
cluster-finding algorithms have based their efficiency mea- 
surements on their ability to detect clusters within artifi- 
cial galaxy data. For example, Postman et al. (1999) and 
Kepner et al. (1999) created artificial galaxy catalogues 
with the same statistical properties as real galaxy cata- 
logues e.g. they matched the surface density of galaxies 
and/or the large-scale clustering properties of the galax- 
ies. They then added artificial clusters to such simulated 
galaxy data. 

We however found this method inadequate since such 
simulated galaxy catalogues cannot fully reproduce the 
hierarchical structure of galaxies in the universe. Specifi- 
cally, the effects of superpositions of structures, variations 
in the field counts, or large-scale structures are hard to 
include in these simulations. Therefore, for each combina- 
tion of RT and LT (see Tables |l| and |), we added a total 
6000 artificial clusters (20 at a time) at random in the 
EDSGC and computed our success rate in detecting these 
artificial clusters above our thresholds (an artificial cluster 
was considered detected if a candidate cluster was found 



by our algorithm within 29c of the original coordinate of 
the artificial cluster; 9c was evaluated at the redshift of 
the artificial cluster). 

In Figure ^, we present our average detection efficien- 
cies as a function of input redshift and richness (using the 
10° X 10° test area discussed in Section |] and in Figure 
m. The error bars are the standard deviation observed 
oetween the different trials of 20 clusters added to the 
EDSGC data at any one time. 

In Figure |^, we show an example of the angular depen- 
dence of our detection efficiency. It is interesting to note 
that our efficiency is strongly correlated with the large- 
scale structure (LSS) in the Universe. This effect is most 
prominent for lower richness clusters, while for the richer 
systems {Rm > 200), it is insignificant (except at the high- 
est redshifts probed by our simulations where the effect re- 
sulted in a loss of ~ 15% of artificial clusters). In addition 
to this interference, our efficiency in detecting low richness, 
high redshift clusters was hindered by the magnitude limit 
of the EDSGC data since the number of potentially visi- 
ble galaxies in these artificial clusters decreases below the 
noise in the background. 

This problem of LSS interference appears to effect both 
over- and under-dense regions of the EDSGC data (when 
compared to the mean surface density of galaxies). For 
example, the most striking example of this LSS interfer- 
ence is seen at coordinates 2.0, 7.5 in Figure ^ where a 
group of clusters - Abell 2730, 2721, 2749, 2755, and 12S 
- has reduce our detection efficiency to almost zero (see 
Figure ^ ). In contrast, we also observe in Figure a low 
detection efficiency near coordinates 8.0, 6.0 which coin- 
cides with a under-dense region in Figure |l|. We believe 
this effect is caused by our assumption of a flat, uniform 
galaxy background as used in the matched filter algorithm 
(see Section 
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In both over- and under-dense regions, 
our assumption of a flat background with the mean surface 
density of the EDSGC is poor and therefore, we suppress 
the overall likelihood of the cluster detection i.e. the model 
of the background around the cluster is not a flat surface 
density oi af — 583775 galaxies per steradian. This ef- 
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feet is further exasperated by systematic plate-to-plate 
uncertainties in the magnitude zero-point of the EDSGC 
photographic plates (Nichol & Collins 1993). 

We present the effective area of the whole EDCCII in 
Table based on our simulation results. These data were 
obtained by summing over the larger EDSGC survey area, 
as defined in Section 0, but weighted by our success rate in 
detecting artificial clusters as computed from the smaller 
test area. The data given in Table ^ illustrate the power 
of this new EDCCII catalogue since we now know the se- 
lection function of an optically-selected cluster catalogue 
to the same accuracy as a X-ray-selected cluster survey 
{e.g. Nichol et al. 1999). This effective area will be used 
below when calculating the space density of clusters (see 
Section |) 

3.2.4. Determining Spurious Detection Rate 

The final use of our simulations was to determine the 
likely spurious detection rate. Again, we have tried to use 
the real galaxy data as much as possible in this analysis 
so as to mimic the real uncertainties in the EDSGC cata- 
logue. This is different from previous attempts to estimate 
the spurious detection rate which have relied on simulated 
galaxy catalogues (see Postman et al. 1996). 

To achieve this goal therefore, we perturbed each galaxy 
in the EDSGC in a random distance (between 29c and 56c 
with a flat distribution) in a random direction from its 
original position. This procedure effectively smoothes the 
galaxy catalogue on these particular scales removing all 
small-scale structure in the catalogue while retaining the 
large-scale features within the catalogue. We then applied 
our matched filter algorithm to these perturbed galaxy cat- 
alogues and calculated the number of clusters that would 
satisfy our selection criteria. This was performed many 
thousands of times to determine the standard deviation in 
our spurious detection rate. 

Our spurious detection rates are shown in Figure ^. The 
number of spurious detections was significant only for high 
redshift Rm — 50 and Rm = 100 clusters. Based on these 
simulations therefore, we restrict ourselves to z < 0.12 for 
Rm < 100 and z < 0.15 for Rm > 100 to ensure that the 
spurious detection rate remains insignificant. For exam- 
ple, for Rm > 200 systems, < 1% of detected clusters are 
likely to be spurious below z = 0.15. We also exclude all 
clusters at z < 0.05 as our thresholds are not accurately 
calibrated below this redshift. In Section ^ therefore, we 
restrict ourselves to Rm > 100 systems in the redshift 
range 0.05 < z < 0.15 and make no correction for spurious 
detections within this richness and redshift range. 

We note here that the results of these simulations are 
in good agreement with empirical determinations of the 
completeness of the EDSGC and EDCC catalogues i.e. 
based on 777 galaxy redshift measurements, Nichol (1993) 
showed that the EDCC was complete out to z ~ 0.13 and 
only contained < 15% contamination from spurious clus- 
ters. At low redshift, Lumsden et al. (1992) also had 
difficulty detecting z < 0.03 clusters in the original ED- 
CCI catalogue because of the large angular size subtended 
by such clusters. 

In addition to randomizing the positions of the EDSGC, 
we also performed simulations which randomly shuffled the 
magnitudes of the galaxies throughout the EDSGC data. 
This resulted in galaxy catalogues with the same statisti- 



cal properties - i.e. same angular clustering and number- 
magnitude relationship - but removed all correlations with 
magnitude. Again, these randomized catalogues were an- 
alyzed with our matched filter algorithm and resulted in 
a very similar result as presented in Figure ^ i.e. the spu- 
rious detection rate was insignificant for lower redshift, 
higher richness systems. We note however, that on average 
we detected fewer rich systems than with the real EDSGC 
data. In other words, magnitude correlations only appear 
to aid in the detection of rich clusters in the data. 

3.2.5. Merging of Catalogues 

Once we have the cluster detections, as a function of red- 
shift and richness, we must then remove duplicate cluster 
detections to produce a final catalogue of unique cluster 
candidates. This was achieved by grouping together all 
systems whose measured centroids were within 2 x 9c of 
each other. We began this process by grouping together 
candidates as a function of their redshift, followed by their 
richness. If duplicates were found, we simply averaged 
their richness and redshift estimates to obtain a final esti- 
mate of the candidate clusters' richness and redshift. 

4. THE SPACE DENSITY OF EDCCII CLUSTERS 

In the section, we present the results of applying the 
match filter algorithm as outlined in Section || to the whole 
EDSGC (as defined in Section ||). However, we must first 
establish a common framework within which to compare 
our results with previous studies. To this end, we will 
quote our results both as a function of Ac/, the cluster 
richness as defined by Postman et al. (1996) for the PDCS 
catalogue, and Rm, which is defined by us. As Rm and Ac; 
are just different richne ss n ormalizations of the match fil- 
ter model (see Section 3.1), it is easy to relate the two 
analytically. This was achieved by integrating a normal- 
ized Schtecter function of richness Rm over the luminosity 



range discussed in Section 3.1 and dividing by L*, the 
characteristic luminosity of the Schtecter function. This is 
equivalent to the original definition of Ac; in Postman et al. 
(1996). Using a = —1.1 (Postman et al. 1996), we show 
that Aci ^ 20.3 is Rm = 100, Ad = 40.6 is Rm = 200 and 
Aci = 81.2 is Rm = 400. We quote here A^; in the PDCS 
V4 passband as this is closest to the EDSGC bj passband. 

In addition to comparing with the PDCS, we wish to 
compare with the original Abell catalogue. We achieve 
this by combining the equality Nji{1.0)/N]i{1.5) ~ 0.7 
(taken from Postman et al. 1996) with Figure 17 of Post- 
man et al. (1996) to obtain an empirical relationship of 
Rabeii ^ 2 X Aci for low rcdshift clusters {z < 0.2). There- 
fore, the three Ac; richnesses given above approximately 
correspond to Abell Richness Classes (RC) 0, 1 and 2 re- 
spectively. These conversions are in good agreement with 
those presented in Lubin & Postman (1996 and references 
therein) and used by others (Holden, private communi- 
cations; Postman, private communications). For the rest 
of this paper, we concentrate on the i?,„ > 100 systems 
(Ac/ > 20 or RC > 0). 

In Tabic ^, we present the cumulative surface densities 
of EDCCII c usters using the methodology outlined in this 
paper. To compute these surface densities, as a function of 
richness, we have summed the number of unique clusters 
with richnesses of > Rm, divided by the effective areas 
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0.05 


0.07 


0.09 
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Table 3 

The effective area (deg^) of the EDCCII as a function of input cluster redshift and richness. These numbers 

WERE EXTRAPOLATED FROM THE EFFICIENCIES WE DERIVED FOR THE SMALLER TEST AREA OF THE EDSGC. 




0123456789 10 
X (Degree) 



Fig. 3. — Our detection efficiency as a function of position in the EDSGC data shown in Figure |]|. Light areas indicate a detection efficiency 
of 100% for our artificial clusters, while dark areas indicate a 0% detection efficiency. Efficiencies near the edge of the area were not measured. 
This plot is for z^st = 0.05 and Rrn = 100 artificial clusters. 




Fig. 4. — Number of spurious detections as a function of redshift for Rm = 50 (top) and Rm = 100 (bottom) clusters in our 10x10 EDSGC 
test area. Number of spurious detections for higher richness clusters were negligible i.e. less than one spurious cluster in the 10°xlO° area. 
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iV(> Rm{Acl),< Zest 

0.12 or 0.15 


) (deg 
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0.4 
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> 100 (> 20.3) 

> 200 (> 40.6) 

> 400 (> 81.2) 


0.41+^-^^ (Ze.t < 0.12) 
0-07At°of2li^est < 0.15) 

0.015 ±0.015 {zest < 0-15) 


3.1 
1.4 
0.015 


4.9 
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Table 4 

Comparison of cumulative surface densities of the EDCCII and the PDCS. The cumulative surface densities of 
the pdcs clusters were calculated using the data in table 4 of postman et al. (1996) using the edccii as a 

ZERO POINT. We excluded all PDCS CLUSTERS WITH a < 3 AND A RADIUS OF > 200 SINCE SUCH CLUSTERS HAVE A HIGHER 
PROBABILITY OF BEING SPURIOUS (SEE POSTMAN ET AL. 1996). We ASSUMED 5.1 deg^ FOR THE SURFACE AREA OF THE PDCS 

SURVEY. Our PDCS surface densities agree with those presented in Figure 21 of Postman et al. (1996). The 
EDCCII data presented in this table has not been corrected for spurious detections. 



Ill Table |. The upper error bars quoted in Table W were 
computed by assuming all clusters detected at a particu- 
lar richness are valid cluster candidates i.e. if we count 
all clusters at that richness regardless of duplicate entries 
(see Section 1^). The lower error bars were computed using 
clusters only detected at that richness i.e. they were not 
detected at any other Rm value. We note that these error 
bars are conservative and should be viewed as boundary 
values. 

In Table |, we present the EDCCII space densities along 
with cluster space density measurements from the PDCS, 
Holden et al. (1999) and the Abell catalogue. For the 
EDCCII, we computed the space densities of clusters, as 
a function of richness(n(i?m)) using the formula 



E 



N{Zest,Rr, 



C v{z) X n{Rm,zesty 



(2) 



which is a sum over all appropriate redshift slices {ZgUce). 
Here, V{z) is the differential cosmological volume (per 



deg2), n{R„ 



is the effective area as a function of 



redshift in Tabled, N{zest,Rm) is the number of clusters 
detected within a redshift slice for a given richness R^ i z\ 
and Z2 are the limits of integration for the redshift slice. 
For Rm < 100, we have summed out to Zest — 0.12, while 
for Rm > 100, we sum to Zest = 0.15. The error bars 
on these measurements were estimated using the same 
methodology as discussed above for the EDCCII cluster 
surface densities. 

5. DISCUSSION 

Recently, Holden et al. (1999) pubhshed 84 redshift 
measurements towards 16 low redshift PDCS clusters. 
From these data, Holden et al. (1999) showed that the 
matched filter redshift estimate for Zest < 0.5 PDCS clus- 
ters had an error of only 5z ~ 0.07; much smaller than 
previously quoted by Postman et al. (1999). Therefore, 
in Table o, we present our estimates of the PDCS space 
density of low redshift clusters (0.2 < z^st < 0.6) using 
the PDCS Zest measurements and the data given in Table 
4 of Postman et al. (1996). We have excluded aU PDCS 
clusters with ct < 3 and a radius of > 200, in the V4 data, 
to minimize the effects of spurious detections (see Post- 
man et al. 1996). Therefore, these space densities may 
be lower than expected since we have potentially excluded 
some real clusters as well. This approach is valid as the 



true error on Zest is now significantly smaller than our red- 
shift slice i. e. a substantial number of clusters will not be 
scattered in, or out, of our sample because of the error in 

^est • 

For the 40 < Ac/ < 80 (RC--1) clusters, we find good 
agreement between the PDCS and EDCCII space densi- 
ties. The same can not be said for the RC~0 systems 
where we differ by over an order of magnitude. We be- 
lieve this discrepancy is the combination of two effects. 
First, the PDCS contains very few < 30 clusters indi- 
cating that the catalogue is possibly incomplete at these 
low richnesses (the PDCS only contains 7 such systems at 
these low richnesses). Secondly, the EDCCII may slightly 
overestimate the space density of low richness clusters as 
the matched filter tends to de-blend rich, nearby systems 
(z < 0.05) into several lower richness clusters. Therefore, 
we will not discuss this inconsistency any further but note 
that this measurement will be important for the next gen- 
eration of large-area CCD cluster surveys since they will 
possess the volume to probe the lower richness systems at 
high redshift. The EDCCII data will provide an impor- 
tant zero-point for such surveys (DeepRange, Zaritsky et 
al. 1997). For RC > 2 systems, the EDCCII appears to 
find a factor of ~ 4 fewer clusters than the PDCS, how- 
ever, the significance of this discrepancy is small given the 
(Poisson) errors on all measurements. 

In addition to checking the matched filter redshift es- 
timates, Holden et al. (1999) also computed the space 
density of PDCS clusters using a new, and completely in- 
dependent, survey selection function than that used by 
Postman et al. (1996) and in this paper. Their results 
are presented in Table ^ and, within the errors, are in 
good agreement with both our PDCS and EDCCII space 
density measurements. Moreover, Holden et al. (private 
communication) have recently re-calculated their measure- 
ments of the PDCS space density, as given in Holden et al. 
(1999) and Table |, but now based on a much larger sam- 
ple of galaxy redshifts (over 700 redshifts in total). They 
now find 18 x 10^^ Mpc"^ for RC=1 systems and 
^ 7 X 10^^ h^"^ Mpc^'^ for RC > 2 systems which is closer 
to the EDCCII space densities presented in this paper. 

In summary, we find good agreement between all three 
of these surveys (PDCS, Holden et al. and the EDCCII) 
which combined span a redshift range of 0.05 < z < 0.6. 
At worst, the difference between the three surveys is i'^^ 
(Table 3). This therefore justifies our original desire to run 
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n(Rm) 



100 < i?™ < 200 (~ 20 < Kci < 40; RC~0) 
200 <Rm< 400 (~ 40 < A^i < 80; RC~1) 
Rm > 400 (~ Aci > 80; RC > 2) 



Cluster Space Densities (10^^ h^^Mpc'^) 
EDCCII PDCS Holden et al. Abell 



83.5^ 



-36.9 

lO.l™ 

9 0+2.5 
'^■'^-2.3 



6.9: 



-3.6 

i2.8t^:^ 

9.411? 



qi O+30.5 



10.4 



+23.4 
-8.4 



11.3 
4.04 
1.58 



Table 5 

Comparison of space density measurements of the EDCCII, PDCS and Abell catalogues. The EDCCII space 

DENSITIES WERE CALCULATED USING EQUATION 1. THE ERROR BARS FOR THE EDCCII CATALOGUE ARE DISCUSSED IN THE 
TEXT, WHILE WE QUOTE 68% ERROR BARS FOR THE HOLDEN ET AL. SPACE DENSITIES AS TAKEN FROM THEIR PAPER. 



the matched fiher algorithm on the low redshift EDSGC 
data since we are now comparing clusters selected in a 
similar way over this entire redshift range. We have thus 
removed one of the main uncertainties associated with the 
PDCS as we do not see a significant difference between 
the low and high redshift cluster populations (as origi- 
nally highlighted by Postman et al. 1996). Moreover, this 
agreement implies that there is little evolution in the space 
density of optical clusters out to z ^ 0.5, in agreement 
with results from X-ray surveys of clusters (see Nichol et 
al. 1997, 1999; Burke et al. 1997; Rosati et al. 1998; 
Vikhlinin et al. 1998; Ebehng et al. 1997, 1999). How- 
ever, we should not overstate this claim, since the error 
bars on all measurements are large. In the future, we will 
need large samples of clusters that span a large range in 
redshift; this should be possible with the next generation 
of cluster catalogues constructed from surveys like DPOSS 
(Gal et al. 1999) and the SDSS (Gunn et al. 1998). We 
also urge the community to adopt one cluster-finding algo- 
rithm so it can be applied to different catalogues (at high 
and low redshift) consistently. 

In Table ^ we present the space density of Abell clusters 
taken from Postman et al. (1996 and references therein). 
The EDCCII space density measurements appear to be 
systematically higher than the Abell catalogue e.g. we 
find ~ 7 times as many RC^O clusters as Abell. However, 
for RC > 0, this discrepancy is much less while the errors 
on these measurements are large. Therefore, we must be 
wary about over-interpretating any claimed discrepancy 
with Abell and simply note that overall, the EDCCII has 
lessened the discrepancy previously claimed to be between 
the high and low redshift cluster populations. 

Our potential disagreement with the Abell catalogue 
could be due to two factors. First, like the EDCCII cat- 
alogue, the effective area of the Abell catalogue could be 
smaller than expected (see Section 3). Secondly, the ED- 
CCII could be finding more clusters than the Abell cata- 
logue at a given richness. The first of these two factors is 
hard to quantify given the subjective nature of the Abell 
catalogue, however, the second factor can be addressed by 
cross-correlating individual clusters in both the EDCCII 
and Abell catalogues. We discuss the latter below. 

In total, we detect 182 of the 324 Abell clusters in the 
EDCCII area, or 56% of them (using a matching radius of 
7.5 arcmins). This is in good agreement with Lumsden et 
al. (1992) who find ^ 70% match-up between their origi- 
nal EDCC clusters and the Abell catalogue. In both cases, 
the percentage of match-ups is independent of richness 



i.e. neither the EDCC or EDCCII appear to have missed 
Abell clusters of a particular richness class (see Figure 8 of 
Lumsden et al.). In addition to comparing richnesses, we 
have also compared the distance estimates of our matched, 
and unmatched, Abell clusters and find little correlation. 
Therefore, the missing 40% of Abell clusters in the ED- 
CCII catalogue appear to be spread evenly over all Rich- 
ness and Distance Classes. Finally, for clusters in common 
between the EDCCII and Abell catalogues, we find no cor- 
relation between the two difference richness estimates i.e. 
Rm and Abell richness or RC. This agrees with Lumsden 
et al. (1992) and Postman et al. (1996) both of whom 
detect a large scatter between their richness estimates and 
the Abell richness estimates. 

In contrast, there are a total of 2109 EDCCII clusters 
detected in the area given in Section |^; a factor of ^ 8 
more than detected in the Abell catalogue over the same 
area (we have excluded the supplementary Abell catalogue 
here and in the above analysis). This discrepancy is less- 
ened when we consider only Rm > 100 systems where we 
find 227 EDCCII clusters (however, we note that the ED- 
CCII only probes to Zest < 0.15, while the Abell catalogue 
contains rich systems out to 2; ^ 0.4). These raw numbers 
reflect the differences in the space densities measurements 
outlined in Table || and highlight that a vast majority of 
the new EDCCII systems are of lower richness (in the ED- 
CCII catalogue). Although, we do note that 143 of the 227 
Rm > 100 systems (63%) are not in the Abell catalogue 
(which corresponds to RC > systems). Again, these find- 
ings agree with the original EDCC cluster catalogue where 
almost 70% of EDCC clusters were new compared to the 
Abell catalogue. These two surveys therefore lend cre- 
dence to the idea (already stated by Abell 1958 and Abell 
et al. 1989) that the Abell catalogue should not be used 
for statistical studies. 
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